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Any time variation in the fine-structure constant alters the ionization history of the universe 
and therefore changes the pattern of cosmic microwave background fluctuations. We calculate the 
changes in the spectrum of these fluctuations as a function of the change in a, and we find that 
these changes are dominated by the change in the redshift of recombination due to the shift in the 
binding energy of hydrogen. We estimate the accuracy with which the next generation of cosmic 
microwave background experiments might constrain any variation in a at z ~ 1000. We find that 
such experiments could potentially be sensitive to \Aa/a\ ~ 10~ 2 — 10~ 3 . 

I. INTRODUCTION 

Physicists have long speculated that the fundamental constants of nature are not constant, but might vary with 
time [Q . Among the possibilities that have received the greatest attention is the time- variation of the fine-structure 
constant a = e 2 /fic. The best laboratory limits on Aa/a give \Aa/a\ < 1.4 x 10 -14 over a period of 140 days 
Limits over a longer timescale can be obtained from astrophysical observations. In particular, spectra from 
high-redshift quasar absorption lines give limits of |Aa/a| < 3 x 10~ 6 at redshifts of z = 0.25 and z = 0.68 ||, and 
| Aa/a | < 3.5 x 10~ 4 for z ~ 3 Q), with a claimed detection at the level of Aa/a = —1.5 ± 0.3 x 10~ 5 for a set of 
redshifts 0.5 < z < 1.6 |]. 

More stringent but also more indirect limits may be placed from geology and cosmology. The Oklo natural nuclear 
reactor yields a constraint of —0.9 x 10~ 7 < Aa/a < 1.2 x 10~ 7 , between a time of 1.8 billion years ago and the 
present raj. Primordial nucleosynthesis gives |Aa/a| < 1.0 x 10~ 4 at a redshift on the order of 10 9 - 10 10 @. 

In this paper, we consider the constraints on Aa/a that could be derived from future observations of cosmic 
microwave background (CMB) anisotropies. Given the plethora of other constraints, is there any reason to examine 
CMB limits on Aa/a? If a is assumed to be constant, then the limits quoted above correspond to \a/a\ < 3.7 x 
10~ 14 /yr (laboratory) 0, \a/a\ < 5 x 10~ 16 /yr (quasar absorption) ||, \a/a\ < 5 - 7 x 10~ 17 /yr (Oklo) |], and 
\a/a\< 1 x 10 _14 /yr (primordial nucleosynthesis) R|. (Here we adopt Hq = 75 km/sec/Mpc, for consistency with 
ref. ||) Our potential CMB limits will not be competitive with any of these. However, in the absence of a particular 
model for changes in a, there is no reason to take a to be constant. Models have been proposed, for example, in 
which a oscillates ||. If the value of a is coupled to a scalar field which evolves on cosmological timescales, then it 
is conceivable that a could vary as a power law in the cosmological scale factor &. One could also imagine models 
in which the scalar field evolves rapidly at early times but later settles into a minimum, producing a fine-structure 
constant which varies at high redshifts, but settles down to a nearly constant value at low redshifts. 

It is useful, therefore, to obtain limits on Aa/a at redshifts z 1. The only limit of this type is provided by 
primordial nucleosynthesis however, that limit is very model-dependent, relying on a particular model for the 
dependence of the neutron-proton mass difference on a. Here we present a much more direct limit, based on changes 
in the spectrum of CMB anisotropies which could be observed by future experiments. 

In the next section, we explain how changes in a alter the recombination scenario, and thus, the CMB spectrum. 
To simplify our discussion, we assume that a has a constant (different) value throughout the recombination epoch; 
i.e., we neglect the possibility that a changes substantially during recombination. In Sec. 3, we calculate the C; 
spectrum for different values of a and explain why our results look the way they do. In Sec. 4, we estimate the limits 
which might be placed on Aa/a at z ~ 1000. We find that the MAP and PLANCK experiments might be able to 
reach sensitivities of |Aa/a| ~ 10~ 2 — 10~ 3 . 
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II. CHANGES IN THE RECOMBINATION SCENARIO 



The fine-structure constant a alters the CMB fluctuations only to the extent that it enters into the expression for 
the differential optical depth f of photons due to Thomson scattering: 



x e n p c<T T , 



(1) 



where <tt is the Thomson scattering cross-section, n p is the number density of electrons (both free and bound) and 

p is the number density of free electrons. The Thomson cross section depends 



x e is the ionization fraction. Thus, 
on a through the relation 



<jt = 87TCt 2 ?i 2 /3?7Jg 



(2) 



The dependence of x e on a is more complicated. Naively, one might expect x e to scale simply with the binding 



energy of hydrogen, which goes as B 



2 /2, suggesting x e (T, a) — x e (T/a 2 ). We will see that this is roughly 



correct, but it is not exact, because the recombination rates depend on a. The reason that x e depends on these rates 
is because it does not track its equilibrium value exactly during recombination. 
Consider the standard ionization equation for Hydrogen JTu]- |11|: 
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where 1Z is the recombination coefficient, (3 is the ionization coefficient, B n is the binding energy of the n th H-atom 
level and n p is the sum of free protons and H-atoms. The Peebles correction factor (C) accounts for the effect of the 
presence of non-thermal Lyman-a resonance photons; it is defined as 



C(a) = 



1 + A 



l + KA(l-x e ) 



1 + A + C l + K(A + 0)(l-x e ) 



(4) 



In the above, K = H^ 1 n p c 3 /87ri/^ 2 (where ^12 is the Lyman-a transition frequency) is related to the ex pan sion time 
scale of the universe, while A is the rate of decay of the 2s excited state to the ground state via 2 photons . Clearly, 
K scales as a~ 6 because ^12 scales as a 2 . Furthermore it can be ascertained that A scales as a 8 jl4j]. To investigate 
/3, one must first use the principle of detailed balance to relate the ionization and recombination coefficients as 



n „ f2irm e kT\ 3/2 

p = n y — p — ) <-*i> 
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while the recombination coefficient can be expressed as 
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where a n i is the ionization cross-section for the (n, £) excited level |I3| . In the above, the asterisk on the summation 
indicates that the sum from n = 2 to 00 needs to be regulated. Physically this comes about due to plasma effects 
which change the ionization and recombination cross-sections (calculated by considering isolated atoms). In essense, 
the summation gets truncated after a certain number of levels. For the present purposes, it suffices to realize that the 
effect of this truncation scheme depends weakly on a and can be neglected [|l5| . 

The a dependence of the cross-section (<r n e) can be summarized as a n i ~ a f(hv/B\), which leads to the equation: 



a \ 



(7) 



This relation is very useful because it allows one to use the temperature parametrizations of 1Z(T) in the literature. In 
particular, TZ(T) can be well fit by a power law of the form T~£. Then from equation ([?]), we see that the a dependence 
of 1Z is just 1Z oc a 2 ( 1+ £). Let the change in a be characterized by A Q = Aa/a ^ 1; then the corresponding fractional 
change in 1Z is 2A Q (1 + £). As it turns out, the results are not sensitive to the precise value of £, which we take to 
be 0.7. Thus, to first order in the change in a, it suffices to consider that 1Z(T) ~ T -0 7 . The ionization equation (||) 
with the change in a can be expressed as 
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dx e 

~~dt 



= c 



Tln v x\ - /3 c ff (1 - x e ) exp 
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where C = (2A Q (1 + £) + l)C(a + Aa) and /3 e ff is the effective ionization coefficient defined as 

f3 eS =pexp(-^(2A a + A 2 a )\ . 



(8) 



(9) 



We have integrated equation (||) using CMBFAST [|T6[ to derive Xq cLS clS <x function of rcdshift for several different 
values of a. The results are displayed in Fig. 1. 
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FIG. 1. The ionization fraction x e as a function of redshift z for the standard scenario (SCDM, Qb = 0.05, h = 0.65) (solid 
curve), an increase of a by 3% (dotted curve), and a decrease of a by 3% (dashed curve). 

The most important feature, the shifting of x e {z) to higher z when a is increased, is easy to understand. Because 
the equilibrium ionization fraction, xf®, is a reasonable approximation to x e , and xf® oc (m e /T) 3 / 2 exp(— B/T), 
which is dominated by the exponential factor near recombination, to a good approximation x e (z) is simply a function 
of z/a 2 (see Fig. 2). 
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FIG. 2. A comparison of the effect on x e of changing a by +3% (dotted curve) with a simple rescaling of the redshift by a 2 
(dot-dashed curve). Solid curve is the original ionization fraction. 

As can be seen in Fig. 2, this scaling is not exact. Two effects spoil it: 1) the factor of (m e /T) 3 / 2 in xf®, and 
2) the fact that x e does not precisely track xf®. Changing a not only changes the energy levels of hydrogen, but 
also all matrix elements and thereby the Thomson cross section and recombination rates. An increase in a increases 
the recombination rates and so equilibrium is more closely tracked. (This can be seen from the fact that the residual 
ionization fraction is smaller for larger a). 

More relevant for the CMB anisotropy is the visibility function, g(z) — e~ T ^dr /dz, which measures the differential 
probability that a photon last scattered at redshift z. The visibility function depends upon x e and <jt through r 
(equation |IJ). The peak of g{z) defines the location of the surface of last scattering and its width determines the 
thickness of the last scattering surface. The finite thickness of the last scattering surface leads to the damping of the 
CMB anisotropy on small scales by smearing out temperature differences on these scales. 

The shape of g(z) is determined largely by x e : around the time of last scattering, the photon mean free path is 
very short until x e — > 0, and the paucity of free electrons makes Thomson scattering rare. Increasing a affects g(z) in 
three ways: first and most importantly, it shifts g(z) to higher redshift because xf® is shifted to higher redshift (by 
the approximate scaling z/a 2 ); second, the larger Thomson cross section increases the opacity by an overall factor, 
which slightly pushes g(z) to lower redshift; and finally, the shape of the g(z) curve is changed because x e more closely 
tracks xf® for larger a. 

Fig. 3 shows the visibility function expressed as a function of conformal time, 77, for different values of a. This is a 
convenient way to display it, because the width corresponds to the comoving damping scale. Increasing a shifts g(z) 
to higher redshift (as explained above); at higher redshift the expansion rate is faster, H(z) oc (1 + z) 3 / 2 , and so the 
temperature and x e decrease more rapidly, making g(z) narrower. The width of the visibility function is predicted to 
scale approximately as l/ot, which is consistent with our results. 
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FIG. 3. The visibility function g(rf) = e -7 dr /drj as a function of conformal time r\ (in Mpc) for the standard scenario (SCDM, 
fib = 0.05, h — 0.65) (solid curve), an increase of a by 3% (dotted curve) and a decrease of a by 3% (dashed curve). The peak 
of g(rf) defines the location of the surface of last scattering and its width defines the thickness of the last scattering surface. 
As can be seen, increasing a moves the last scattering surface to higher redshift (smaller conformal time) and decreases its 
thickness. 

Are there any other potential effects on the CMB due to a variation in a? One completely negligible effect is the 
change in the He recombination scenario due to the change in the binding energies of He atomic levels. Another effect 
is the change in the variation of the matter temperature with time. Specifically, the matter temperature variation 
consists of adiabatic cooling due to the expansion of the universe and the cooling due to Thomson scattering. The 
change in ut changes the latter. However, the matter temperature accurately tracks the radiation temperature until 
very late (1% difference at z ~ 500) and hence this effect has no consequences for the present purposes. 

III. CHANGES IN THE CMB FLUCTUATION SPECTRUM 

We have integrated the changes in the differential optical depth due to a variation in a into CMB FAST |l6|. The 
results are shown in Fig. 4 for a ±3% change in a. Two separate effects may be noted from the results. One, for 
an increase in a, the peak positions in the spectrum shift to higher values of I. Two, increasing a causes the values 
of Ci to systematically increase. Conversely, a decrease in a shifts the peaks to lower values of I and decreases their 
amplitude. 
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FIG. 4. The spectrum of CMB fluctuations for the standard scenario (SCDM, fij, = 0.05, h = 0.65) (solid curve), an increase 
of a by 3% (dotted curve), and a decrease of a by 3% (dashed curve) 

To understand the first feature, a qualitative understanding of the position of the peaks is necessary. Using £ p 
to denote the position of a peak, rg(z) for the angular diameter distance and r s {z) for the sound horizon, one can 
write jnj £ p ~ rg(zi s )/r s (zi s ), where zi s is the redshift of the surface of last scattering. Increasing a increases the 
redshift of the last scattering surface, as seen in Fig. 3. A higher redshift at the last scattering surface corresponds 
to a smaller sound horizon and thus, a higher value of I. Decreasing a has the opposite effect: the redshift of last 
scattering decreases, producing a larger sound horizon at last scattering, and thus a smaller value of I for the peaks. 

The increase in the amplitude of the peaks with increasing a derives from two separate effects. The amplitude 
of the first peak is quite sensitive to the magnitude of the integrated Sachs- Wolfe (ISW) effect. If a mode enters 
the horizon when radiation still makes a significant contribution to the energy density, the decay of the gravitational 



potential leads to the blueshift of photons |18|. This effect has been dubbed the "early ISW effect" to distinguish 
it from the decay of the gravitational potential at late times in models which become dominated by curvature or a 
cosmological constant. An increase in a pushes recombination to a higher redshift, resulting in a larger early ISW 
effect and, thus, a larger amplitude of the first peak. The early ISW effect is felt most strongly around the scale of 
the sound horizon at last scattering. For the SCDM model we have considered, this is around 100 Mpc or £ ~ 100. 
By I ~ 500, the effect of early ISW contributions is negligible. 

Beyond the first peak a second effect is dominant: diffusion damping of CMB fluctuations due to the finite thickness 
of the last scattering surface (see Fig. 3). Because the last-scattering surface is not infinitely thin, the anisotropies 
seen today are an average over a region of finite thickness defined by the visibility function. This leads to damping of 
small-scale anisotropies, given by a photon diffusion damping factor averaged over the visibility function [ fL9| , 

/>oc 

D(X) = dzg(z)exp(-X 2 D (z)/X 2 )^eM-^D(zis)/X 2 )- (10) 



The characteristic damping length Xd is set by the width of the visibility function. (The multipole damping scale is 
given approximately by Id ~ /Xd)- As explained earlier and shown in Fig. 3, the comoving damping length 

decreases with increasing a. Thus, an increase in a decreases the effect of damping, and the power spectrum at large 
I increases with increasing a, as seen in Fig. 4. 
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IV. LIMITS ON VARIATIONS IN THE FINE-STRUCTURE CONSTANT 



From the analysis presented in sections [0| and [II, it is clear that a variation in a has a substantial effect on the 
CMB fluctuation spectrum. The aim of this section is to obtain a quantitative measure of the limits put on a by an 
ideal CMB anisotropy experiment. This can be accomplished through an analysis of the Fisher information matrix. 
If our estimate of the cosmological parameters is very close to the true values, then the likelihood function (£) 
can be expanded about its maximum as 

£~£ m exp(-iV^), (11) 
where F^ is the Fisher information matrix, defined as [po[ 

1 fdCA fdCA , . 



! 2 ACf \dO*J \d0j 



In equation (|T2|), the quantity ACi is the error in the measurement of Ci. From the Gaussian form of C, the covariance 
matrix is seen to be F -1 . In particular, one can define the standard deviation for each parameter 9i as of = (F~ 1 )a. 
The cosmological parameters (9i) that need to be determined from the measured fluctuation spectrum are taken 
to be the Hubble parameter (h), the number density of baryons (parametrized as f^ft 2 ), the cosmological constant 
(parametrized as Q.{Ji 2 ), the effective number of relativistic neutrino species (N v ), the primordial helium mass fraction 
(Y p ), and the fine-structure constant (a). We make the assumption that the experiments are limited only by the cosmic 
variance up to a maximum £, denoted by £ rn ax- This assumption is an oversimplification, but it provides a rough 
upper bound on the possible limits on Aa/a from future CMB experiments. 

The fiducial models used for the present work are a standard cold dark matter model (SCDM) and a CDM model 
with a cosmological constant (ACDM). Both models have fl b h 2 = 0.02, h = 0.65, Y p = 0.246, and N u = 3.04. (Note 
that various higher-order effects, most notably the slight heating of the vv pairs by electron-positron annihilation, 
increase the effective value of N v to 3.04 from its canonical value of 3 |pl[ ). In the ACDM model, Q\ is taken to be 
0.7. We use an adiabatic, scale invariant initial power spectrum and constrain the cosmology to be flat in keeping 
with the standard inflationary paradigm. For each of these two models, we consider two limiting cases regarding prior 
constraints on the unknown parameters: first, no prior constraints at all, and second, a "best-case" set of limits on 
the unknown parameters using priors | p2| . In the latter case, we take, as 1 — a limits, h = 0.65 ± 0.05 from current 
observations, and fif,ft, 2 = 0.02 ± 0.002 and Y p = 0.246 ± 0.001 from Big-Bang nucleosynthesis j23|. For this case, we 
also fixed N v to be exactly equal to 3.04. 

The required derivatives of the C/s were calculated by two-sided finite differencing for each parameter, while the 
rest were kept fixed. We verified that the changes in the results obtained were less than 10% when the variation of 
the parameters was halved. The results are shown in Fig. ^|in terms of the ratio a a /a, where a a is the 1 — a accuracy 
measure obtained from the Fisher matrix analysis. We see that the estimated upper limits on |Aa/a| vary from about 
10~2 for £ max ~ 500 - 1000 down to 10" 3 for £ mnr > 1500. 
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FIG. 5. The estimated accuracy with which a can be constrained by a cosmic variance limited CMB anisotropy experiment, 
as a function of the maximum angular resolution given by £ max . The dotted curve is the result of including priors as explained 
in the text, while the dashed curve is for the case without priors. 

These results suggest that future CMB experiments (MAP and PLANCK) might be able to constrain any variation 
in the fine-structure constant to less than 1CP 2 — 10~ 3 . This is a weaker constraint than can be obtained from current 
quasar absorption studies, but the CMB limit would apply at a much higher redshift (z ~ 1000). It represents a much 
more direct and reliable constraint than the only other limit at z ^> 1, available from Big Bang nucleosynthesis 0. 
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